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Abstract 

Neutron stars are one of the most exotic objects in the universe and a unique 
laboratory to study the nuclear matter above the nuclear saturation density. In 
this work, we study the equation of state of the nuclear matter within a relativistic 
model subjected to a strong magnetic field. We then apply this EoS to study and 
describe some of the physical characteristics of neutron star, especially the mass- 
radius relation and chemical compositions. To study the influence of a the magnetic 
field and the hyperons in the stellar interior, we consider altogether six solutions: 
three different values of magnetic field to obtain a weak, a moderate and a strong 
influence, and two configurations: a family of neutron stars formed only by protons, 
electrons and neutrons and a family formed by protons, electrons, neutrons, muons 
and hyperons. In all cases the particles that constitute the neutron star are in (3 
equilibrium and zero total net charge. In the end, we see that although a magnetar 
can reaches 2.50A/q, arises a natural explanation of why we do not know pulsars 
with masses above 2.0Mq. 
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1 Introduction 

Neutron stars are compact objects maintained by the equilibrium of gravity and the 
degenerescence pressure of the fermions together with a strong nuclear repulsion 
force due to the high density reached in their interior. Since we do not know yet 
the precise and detailed structure and composition of the inner core of a neutron 
star, many models have been used to describe it. In the literature we can find 
some standard ones: hadronic neutron stars, quark stars, strange stars and hybrid 
stars [TH3]. 

In the present work we study a hadronic neutron star constituted by nucleons 
and hyperons and subject to a strong magnetic field. The presence of hyperons 
is justifiable since the constituents of neutron stars are fermions. So, according to 
the Pauli Principle, as the baryon density increases, so do the Fermi momentum 
and the Fermi energy. Ultimately the Fermi energy exceeds the masses of the 
heavier baryons pQ . On the other hand, some strange objects like the soft gamma- 
ray repeaters and anomalous X-ray pulsars can be explained assuming that these 
objects are neutron stars subject to a strong magnetic fields on their surface. These 
objects are called magnetars [4 . Although the magnetic field of the magnetars 
do not exceed 10 15 G in their surface, it is well-accepted in the literature that 
the magnetic field in the core of the neutron stars can reach values greater than 
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10 18 G [5] [ 6 ] . Due to the large densities in the neutron star interior, we do not expect 
any significant influence of the magnetic field till it reaches values of the order of 
10 18 G. 

This paper is organized as follows: we make a review of the formalism of the 
non-linear Walecka model (NLWM) in the presence of a magnetic field. Then we 
present the numerical results showing how the presence of hyperons and a strong 
magnetic field affects the EoS and the chemical composition. Finally we study how 
these terms alter the macroscopic mass-radius relation of the neutron stars and 
compare our results with those found in literature. 

2 The Formalism 

The total Lagrangian is given by [TIE] : 

c = ^2c b + c m + ^2ci + c B) (1) 

b I 

where b stands for the baryons, m for the mesons, I for the leptons, and B for the 
electromagnetic field itself. The sum in b can run over the eight lighter baryons and 
in I over the two lighter leptons. Explicitly, in the presence of a electromagnetic 
field the Lagrangian is: 

C b = ^bhuiid" - eA» - g v<b u» - g P ,bhbP») - (M b - g s , b (r)}^ b , (2) 

+ \™%P^) - \P^ V - t^ 3 - \\a\ (3) 
h = Ulu(id» -eA^-mi}^, (4) 



1 
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where ^S b and tpi are the baryon and lepton Dirac fields, respectively. The baryon 
mass and isospin projection are denoted by M b and /3b respectively. The mass of 
leptons are mi and the electric charge of the particles is given by e. The antisym- 
metric mesonic and electromagnetic field strength tensors are given by their usual 
expressions: = d u Lo v - P pu = d p p v - d v p p and = d p A v - d v A ll . 

The 7^ are the Dirac matrices [9]. The strong interaction couplings are denoted 
by g, and the meson masses by to, all with appropriate subscripts. The second 
subscript at the g constant is due to the distinctive coupling of hyperons with the 
mesons. In this work we assume that g Sy H = 0-7g s ,N', 9v,H = 0.783^^ an d g P ,H — 
0.783g Pi 7v |10j . where H denotes hyperons and N nucleons. The hadronic part of the 
Lagrangian is the called NLWM. The leptons are included in the total Lagrangian 
density as a non-interacting Fermi gas in order to account the /? equilibrium in the 
star. 

To solve the equation of motion, we use the mean field approximation, where 
the meson fields are replaced by their expectation values, ie: a — > (cr) = Co, ^ ^ 
<V (w M > = ujq and -> S 0p (p M > = p - 
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In this work we use a GM1 parametrization which can describe the most 
important properties of nuclear matter and reproduce the macroscopic properties 
of the neutron stars consistent with those observed in nature. The GM1 parameters 
are: 



Set 


{g s /m s f 


(g v /m v y 


{g P /m p f 


k/M n 


A 


GM1 


11.785 fm z 


7.148 fm z 


4.410 fm? 


0.005894 


-0.006426 



Table I: Values of GM1 parametrization. 



This parametrization is fixed so that the incompressibility of nuclear matter 
K = 300 MeV, and the nuclear saturation density no = 0.153/to -3 . The masses 
of the baryon octet are M N = 939 MeV (nucleons), M A = 1116 MeV, M E =1193 
MeV and M s = 1318 MeV. The meson masses are m s = 400 MeV, m v = 783 
MeV and m p = 770 MeV. And the masses of the leptons are m e = 0.511 MeV and 

= 105.66 MeV. Applying the Euler-Lagrange in equation Jl]) in the absence 
of an electric field, the equation of motion in the mean field approximation for an 
arbitrary baryon becomes: 

bb(*0° - 9v,iW> - 9 P ,bhbPo) ~ lj {id 3 - eA>) - M b *]* = 0, (6) 

where 



M; = M b - 9s>b a , (7) 

is the baryon effective mass. 

For an uncharged particle eA^ is always zero. The quantization rules are: id 
= E and id 3 = where k 3 is the momentum in j direction. Setting E — g v ^o ~ 
9p,bl3bPo — E* , we have the following equation of motion written in a block matrix: 



0. (8) 



(E* - M b *) <j.k \ / u A 

a.k - {E*p + M*) J \ u B 

This is an eigenvalue equation, which can be solved as the free Dirac equation 
for an effective mass and energy, whose solution is: 



E = \Jk 2 + M* 2 + g Vtb oj + g p ,bhbPo = M, (9) 



where fj, is the chemical potential. 

For a charged baryon, the Dirac equation assumes the following form: 



(E*-M*) -(T.(k-eA) \ f u A 
a.(k - eA) - (E* - M* b ) \u B 



0. (10) 



To produce a constant magnetic field in the z direction we choose: A2 = A3 = 0, 
A\ = —By. The solution of this eigenvalue equation is well-know in the literature 
EH2]: 



E* = JM b * 2 + fcf + 2u\e\B, 



E = X /M* 2 + k 2 + 2v\e\B + g Vjb cj + g P , b hbpo = P, (11) 
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where the discrete parameter v is called Landau level (LL) and fi is the chemical 
potential. The first LL, v = 0, is non degenerate and all the others are two-fold 
degenerate. For the leptons, since they don't feel the strong force: 



Ei = Jmf + k% +2v\e\B = tn. (12) 
And the expected values for the mesons are: 
u> 



y ^A n ub + y 9_v± nbc (13) 

ub cb 



\^ 9s,b ub sr^ 9s,b he 1 K 2 1 A a h a \ 



ml s ^— ' m„ s 2 m„ 6 m„ 

u6 cb * s s 



Po = Y.H nUb hb+Y.H nCbI ^ ( 15 ) 

where n cb and n" h are the number density of the "charged baryons" and "uncharged 
baryons" respectively [T31III], and nf 1 and n" b are called scalar density for the 
charged and uncharged baryons [7] . In T = Qj they are given by: 

eli? 
{2ny 



\ ' w J —K} v 



n 



* = v „ M / fe/ M ^ fc * d9) 



The summation in v in the above expressions ends at v maxi the largest value of 
v for which the square of Fermi momenta of the particle is still positive and which 
corresponds to the closest integer from below defined by the ratio: 

2 _ M *2 

v max < i ip i charged baryons (20) 
2\e\B 

2 2 

/i — mf 

v m ax < 2 | e |^ ■ leptons (21) 
Now we couple the equations imposing ft equilibrium and zero total net charge: 

j"6i = A«n - GifXe, Me = M/ii /J & b^b + e W = °, ( 22 ) 



1 This can justified since the Fermi temperature of the neutron stars is too high com- 
pared to its own temperature [15] . 
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where fi^ and are the chemical potential and electric charge of the i-th baryon, 
and n n , fi e and are the chemical potential of the neutron, electron and muon 
respectively. 

The energy density of the neutron star is: 

ub cb I m 

where the energy densities for the uncharged baryons, charged baryons, leptons and 
mesons have the following forms: 



e„6 = -5 / \/M* 2 + ^ 2fc2 *» ( 24 ) 



e «* = E jf ' >/M 6 * 2 + ^ + 2i/|e|Bdfe a 
E / ' \/mf + k 2 z +2v\e\Bdk z 



(25) 



e m = ^^crg + im^wg + ^m 2 pP l + ^KO% + ^| Acr o- ( 27 ) 
To find the pressure, we use the second law of thermodynamics, that gives: 

B 2 

V = E ^ lU ' 1 ~ e + g^' ( 28 ) 

i 

where the sum runs over all fermions. Note that the contribution from electromag- 
netic fields should be taken into account in the calculation of the energy density 
and the pressure. 



2.1 TOV equations and the density-dependent magnetic field 

The magnetic field of the surface of the magnetars are of order of 10 15 G, but 
can reach more than 10 18 G in their cores. To reproduce this behaviour we use a 
density-dependent magnetic field given by [SJ[T21[TI>] : 



B(n) = B suf + B 



1 — exp 



(29) 



gsuf j s ma g ne tic field on the surface of the neutron stars, taken as 10 15 G, n 
is the total number density, n = ^n b , Bq is the constant magnetic field. The 
two free parameters /3 and a are chosen to reproduce a weak magnetic field bellow 
the nuclear saturation density, and a quickly growing when n > no, in such a way 
that B(n) ~ Bq when n > 6no- To reproduce this behaviour, we have set /3 — 
—5.76 • 10~ 3 and a = 3. Now B is replaced by B(n) in the term B 2 /Stt in our EoS. 
To finish our analytical analysis, we write the TOV |17j equations: 



(30) 
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dp 
dr 



Ge(r)M(r) 



1 



p(r) 



e(r) 



1 



Aitp(r)i 



M(r) 



1 



2GM(r) 



(31) 



which are the differential equations for the structure of a static, spherically symmet- 
ric, relativistic star in hydrostatic equilibrium. The equation of states developed 
in this work are used as input for these equations. Since the distribution of the 
magnetic field is not spherical, we are aware that the TOV equations can only be 
used as an approximation in the present study. 

3 Results and discussion 

We consider two families of neutron stars: one containing just protons, electrons 
and neutrons, which we call "Atomic Stars" denoted by the letter A in the legends, 
and other containing protons, electrons, neutrons, muons and hyperons, which we 
call "Hyperonic Stars" denoted by the letter H in the legends. In the results we 
also include the crust of neutron star through the BPS EoS [IB], but always taken 
into account the contribution of the magnetic field through the term B(n) 2 /8ir in 
the EoS. 
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Fig. 1: EoS for two atomic and two hyperonic stars obtained with different 
values of the magnetic field. The straight line corresponds to the 
causal limit, for which e = p. 



We choose three values for the magnetic field: 2.2 • 10 17 G, 3.1 • 10 18 G and 6.6 • 
10 18 G, to produce a weak, a moderate and a strong influence. We also include here 
some theoretical and observational constrains. First, all our EoS are causal and 
obey the Le Chantelier principle, ie, the quantity dp/de lies between and 1. We 
plot the numerical results of four EoS in fig. 1. 

As we can see from fig. 1, the presence of hyperons softens the EoS more 
than the influence of the magnetic field can stiffen it. No matter how strong is the 
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magnetic field in the interior of the magnetar, the EoS of an atomic star is always 
stiffener than a hyperonic one. We also can see that all our EoS are causal (the "c" 
line is the causality limit). 

From figs 1 and 3 we note that both the EoS and the fraction of particles 
Yi = rii/n are not affected significantly by a magnetic field about 2.2 • 10 17 G. The 
reason is that a field of this magnitude is too weak to contribute to the final EoS 
and to the fraction of particles. We also see that the magnetic field affects more the 
hyperonic stars than the atomic ones. As the hyperonic stars are softer than the 
atomic ones, they are therefore more sensitive to the presence of the magnetic field. 
Also the hyperonic stars are denser than the atomic stars, with a bigger central 
density n c . Moreover, for a fixed value of density n, the pressure p of the hyperonic 
stars are smaller than the atomic ones. So, as the magnetic field couples to the 
number density through the equation (|29[) . the contribution of the magnetic field 
is always greater in hyperonic than in atomic stars. We show this result plotting 
the pressure p in function of n in fig. 2. 

4 I 1 T« 1 1 f ' 1 

A, B-6.6. 10 ®G 

A, B = 2.2 . 10 'G 

H, B = 6.6. 10 °G 
H, B-2.2 . 10 G 

3 - 



"e 2 - 

Oh 




1 f- 1 1 1 1 

0.2 0.4 0.6 0.8 1 

n (fm~ 3 ) 



Fig. 2: Pressure as function of number density for four stars. 

Figures 3, 4 and 5 show how the fraction of particles change with the increase of 
the magnetic field. We can see that for a considerable strength field, the appearance 
of charged particles is favoured at low densities due to their dependence on the 
magnetic field, as expected from eq. (|T7|) . Another change is the strange behaviour 
of the particles in the presence of a considerable magnetic field. For weak fields, the 
population of a kind of particle is always well-behaved, while for the strong ones 
many kinks appear. 
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The reason is that in the absence of a magnetic field the number density of a 
determined kind of particle grows smoothly with the momentum. When a magnetic 
field is present, there is also a dependence of the discrete LL. For a weak magnetic 
field a lot of Landau levels are occupied, but for a strong magnetic field, just a few 
of them arc filled. So the orbit normal to the z direction is tightly quantized. This 
effect is more evident in the hyperonic star. 

Fraction of particles Yi as a function of number density n. 
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Fig. 3: (a) atomic stars and (b) hyperonic stars for a magnetic field of 2.2 . 
10 17 G 




Fig. 4: (a) atomic stars and (6) hyperonic stars for a magnetic field of 3.1 
10 18 G 
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Fig. 5: (a) atomic stars and (6) hyperonic stars for a magnetic field of 6.6 . 
10 18 G 

Each nozzle in a slope of a determined particle indicates that the density is 
high enough to create another Landau level. For a high density there are so many 
Landau levels available that the distribution approaches to the continuous, while 
for a weak magnetic field, even in a low density there are several LL, so there is no 
significant difference with the case in the absence of the magnetic field. One can 
notice that for hyperonic stars, at densities of the order of 0.8/m -3 the neutron 
is no longer the most important constituent. From this point, the A hyperon 
dominates in the region of high densities. 

It is hard to compare quantitatively our results with other ones existing in the 
literature due to the many different parametrizations. In the absence of hyperons, 
we see that our EoS is stiffener than those presented in [7J[5] since the authors of 
those references do not consider a density-dependent magnetic field. With relation 
to the fraction of particles Yi , the fraction of protons in a low density region is more 
favoured in [7] than in our work. When hyperons are present, for a magnetic field 
of 2.2 • 10 17 G there is no difference in both, the EoS and in Yi when compared with 
the results for zero magnetic field pQ. When a strong magnetic field is applied, we 
see that our parametrization do not prevent the hyperon formation in a low density 
region as showed in |19j . Also, in our work, the influence of Landau quantization 
is much more evident even for a smaller value of the magnetic field. 

In order to validate our EoS, we have to solve the TOV equations and check 
if the results agree with observational constraints. The star masses cannot exceed 
the maximum theoretical neutron star mass of 3.2 solar masses |20j . The EoS 
has to be able to predict the 1.97 solar masses neutron star [21] and to be in 
agreement with the redshift measurements (z) of two neutron stars. A redshift of 
z = 0.35 has been obtained from three different transitions of the spectra of the 
X-ray binary EXO0748-676 [52]. This redshift corresponds to M/R=0.15 M Q /km. 
Another constraint on the mass-radius ratio comes from the observation of two 
absorption features in the source spectrum of the IE 1207.4-5209 neutron star, 
with redshift from z = 0.12 to z = 0.23, which gives M/R=0.069 M Q /km to 
M/R=0.115 MQ/km [23] . Besides these constrains, the stars with central density 
above that of the maximum mass stars are mechanically unstable pQ. Due this fact, 
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the S° hyperon is not present in the neutron star interior (the density required to 
create it is too higlQ). 

Solving the TOV equations for the EoS we obtain the results presented in Table 

2 . 





Mass (Mq) 


Radius (km) 


n c (/m~ 3 ) 


B (xl0 i8 G) 


Max. Mass (A) 


2.50 


12.30 


0.701 


6.6 


Max. Mass (A) 


2.42 


12.15 


0.798 


3.1 


Max. Mass (A) 


2.39 


12.10 


0.840 


0.22 


Max. Mass (H) 


2.28 


11.82 


0.796 


6.6 


Max. Mass (H) 


2.09 


11.69 


0.907 


3.1 


Max. Mass (H) 


2.01 


11.84 


0.952 


0.22 



Table II Neutron stars properties computed from the six EoS used as input to the 
TOV equation. From the central density, we see that the hyperonic stars are 
denser than the atomic ones. 



The fact that the EoS of hyperonic stars are always softer than the atomic ones 
reflects in the maximum mass of these stars. Further, from Table II we can see that 
the hyperonic stars are denser than the atomic ones, as said before. Ultimately, 
there is a curious behaviour of the radius of the maximum mass hyperonic stars. 
While in the atomic stars the radius grow with the magnetic field, in the hyperonic 
ones the radius has no visual dependence with it. We plot the TOV solutions in 
fig. 6 to compare our results: 

3 



2.5 
2 

O 
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s 

1 
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7 8 9 10 11 12 13 14 15 16 

R (Km) 

Fig. 6: Mass-radius relation for two atomic and two hyperonic stars with 
different values of magnetic field. The straight lines are the observa- 
tional constrains 

2 indeed the H° hyperon appears in an insignificant quantity Y s a = 10 -5 at 1.2/m -3 . 
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In fig. 6 we can see that all our models are in agreement with the previously 
discussed constrains [20-23 . The inclined straight lines are the constrains of the 
measured redshift while the horizontal one is the 1.97 Mq pulsar. Every single dot 
in the curves is a possible neutron star. We can see also that a magnetar can reach 
a mass of 2.50Mq. However, as we do not known neither expected to discover 
any white dwarf with mass above IAMq due to the Chandrasekhar limit 0, the 
same behaviour appears to the neutron stars with mass above 2.0Mq due to the 
Oppenheimer-Volkoff limit if we believe that hyperons and muons existing on their 
interior. The maximum possible max is very close of the 1.97 Mq known neutron 
star. The discovery of a neutron star with mass above 2.0Mq will imply that 
hyperons are absent or the star is a rare super magnetic neutron stars, as shown in 
fig- 6. 

4 Conclusion 

In this work we consider a hadronic neutron star compost by hyperons submitted 
to a strong magnetic field. We see that while the presence of hyperons reduce the 
maximum mass by softening the equation the state (EoS) [Tl l2"4"H2"7] , the presence of 
a density-dependent magnetic field tends to increase the maximum mass stiffening 
the EoS [12j|T6] ,Our study shows also that although the influence of a strong mag- 
netic field stiffener the EoS, a hyperonic magnetar still has a softer EoS compared 
with a common atomic stars. We also see how the magnetic field can change the 
chemical composition of neutron stars due the Landau quantization and offer an 
explanation about the non-existence of neutron stars with mass above 1.97 Mq. 
We conclude this paper by indicating that the effect of a strong magnetic field has 
to be taken into account in a description of magnetars, mainly if we believe that 
there are hyperons in their interior. In this case, the influence of the magnetic field 
can increase the mass by almost 15%. 
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